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ABSTRACT 

M 82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In 
addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford 
better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution 
and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more 
extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others 
with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive 
spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical 
conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, 
and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind. 
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1. Introduction 

M 82 is the closest specimen of sta rbursts with superwin ds, at 
a distance of only (3.9 + 0.3) Mpc (ISakai & Madordll999l) . Its 
almost edge-on orientation is ideal for the study of the super- 
wind, and more generally the interaction of galaxies with the 
intergalactic medium and the role of the gas and dust dynamics 
in the self -regulation of starbursts. 

The exceptional star formation activity of M 82 is thought to 
be caused by its tidal interaction with the M 8 1 galaxy group, 
redistributing the gas within the disk of M 82 toward the cen- 
tral regions where it accumulates, and expelling atomic gas in 
giant streams c onnecting the disk to the intragalactic medium 
(lYun et aO 19941) . M 82 is an evolved starburst, in the sense that 
its activity has been sustained (episodically) over much of its 
lifetime. Coupled stellar population synthesis and photoioniza- 
tion models reveal two major episod es of s tar formation over 
the past 10 Myr (Forster Schreiber et al. 2003). Numerous young 
supernova remnants have been detected ( Muxl ow et al.l H994I) . 
Another peak of star formation about 1 50 Myr ago is inferred 
from the age-dating of stellar clusters dKonstantopoulos et al.l 
2009), and 200 super star clusters, signposts of starbursts, have 
been found in M82 ( Melo et al.l 120051) . spanning a wide age 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led consortia, and important participation of NASA. 



interval from 60 My r to several Gyr dGallagher & Smithl[J999t 
Ide Griis et al.ll2001l) . 

IWard et all (l2003h found that more than half the total molec- 
ular gas o f the central regions is warmer than 50 K, and 
iMao et alj d2000t) found that CO transitions up to J - 7 - 
6 trace diffus e gas abundantly illuminated by UV radiation. 
iPanuzzo et al.l (2010) modeled CO transitions from 7 = 4-3 
to / = 13-12 in the central 43", observed with the SPIRE 
FTS spectrometer, and they suggest that turbulent heating caused 
by stellar winds and supernovae may be the dominant exci- 
tation process for the highest transitions, tracing gas at about 
500 K. The molecular gas properties thus reinforce the notion 
that M 82 is an evolved starburst, which has been active for long 
enough to completely disrupt the cold component of the inter- 
stellar medium. 

The mechanical energy simultaneously released by a 
large collection of massive stars and supernovae is able 
to drive a superwind detected in both X-rays emitted by 
shocked gas (Strickla nd etal.1 [1997) and optical lines ex- 
cited b y a combination of photoionizatio n and shock ion- 
ization (Shopbell & Bland-Hawthorn 1998). Optical emission 
around the wind al so partly issues from a reflection nebula 
(Scar rott et alJ ll991). indicating the presence of dust. Abundant 
molecular gas is associated with both the superwind (although 
probably not spatially coincident but rather outside t he bicylin- 
drical outflow) and tidal streams dWalter et al.ll2002l) . and PAH 



Roussel et al.: SPIRE imaging of M 82: cool dust in the wind and tidal streams 



emiss ion is found in a very extended halo (Engelbrach t et al.l 

This paper studies the spatially-resolved properties of the 
cool dust detected by SPIRE. Since PACS photometric observa- 
tions had not been scheduled at the time of writing, and Spitzer 
MIPS images are saturated, we sample the far-infrared spectral 
energy distribution only longward of 250 /mi, where it can be 
safely assumed that the emission is dominated by grains in ther- 
mal equilibrium. SCUBA maps have higher angular resolutions 
than SPIRE maps; their sensitivities are however too low to reli- 
ably detect any emission outside the disk of M 82, even when all 
archival data are combined ( Leeuw & Robson 2009). In contrast, 
bright filaments and diffuse emission are detected out to very 
large distances with SPIRE. We thus focus on the dust content 
and physical properties outside the main body of M 82, where 
SPIRE maps provide truly unique information. 



2. SPIRE observations, data reduction, and analysis 

M82 was observed with the SPIRE photometer arrays onboard 
Herschel (Pilbratt et al. 2010) during the science demonstration 
phase. The SPIRE instrument, its in-or bit perform ance, and sci- 
entific capabilities are described by Griffin et al. (2010). Its re- 
solving power enables a linear resolution of 350 to 700 pc in 
M 82, from 250 to 500 /mi (1 kpc corresponds to 53"). 

SPIRE maps of 22' x 22' were acquired in nominal scan 
mode, with a total of 4 scans. The d ata were processed with HIPE 
3.0 up to level 1, as described by iPohlen et alj d2010h . except 
without the thermal drift correction. The data were then input 
to the Scanamorphos software (Roussel, in preparation), which 
produces maps after subtracting the thermal drift and the low- 
frequency noise of each bolometer, purely by using the available 
redundancy. The resulting maps have ~ 30% smaller sky stan- 
dard deviations and ~ 3% larger beam areas. 

We also derived SPIRE point response functions (PRF) Q 
from a nominal 8-scan observation of Uranus, with a field of 
view of 20' x 20', that we processed in the same way and with the 
same orientation along the first scarfj. This observation is deeper 
than that of Neptune - made in bright-source mode - used for the 
official PRF release, and thus allows detection of the secondary 
lobes and diffraction spikes. Uranus is slightly more extended 
than Neptune (3.53" at the time of observation^, as opposed to 
2.26"), but suitable for sources larger than 70pc in M82. To 
create the PRF, all the emission outside the beam pattern was 
masked out, and the few remaining background sources were 
approximated by Gaussians and subtracted. 

All maps have a pixel size of 4.5" at 250/mi, 6.25" at 
350/mi, and 9" at 500 //m, i.e. approximately one fourth of the 
beam full width at half maximum (FWHM). 

Photometry: Since the flux calibration of SPIRE is strictly valid 
for point sources only, we used our empirical beams (corrected 
for the finite size of Uranus) to estimate beam areas of 468, 
850, and 1763 square arcseconds, respectively at 250, 350, and 
500 /mi, and obtain extended-source flux calibration. They are 
larger than the current official estimates (based on Gaussian fits, 
i.e. not including sidelobes and spikes) by 10 to 20%. 



We measure global flux densities of (457. ± 2.) Jy at 250/mi, 
(155. ±2.) Jy at 350/mi, and (49.6 + 0.9) Jy at 500 /mi (error bars 
not including systematic flux calibration uncertainties of 15%, 
nor beam size uncertainties). Th ey were multiplied by the rec 
ommended corrective factors (see Pohl en et 



implied b 
aIl l20Toh . 



1 This term designates the result of sampling the point spread func- 
tion (PSF) of the optics by the combination of the observation pattern 
and the map pixel grid: it is an empirical realization of the PSE 

2 files available at http://www2.iap.fr/users/rousseI/spirebeams/ 

3 using the HORIZONS ephemeris service 
http://ssd.jpl. nasa. gov/horizons. cgi 
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Source decomposition: The central regions of M82 have a 
complex structure. They contain multiple bright infrared sources 
that can be most easil y separated in m id-infrared images 
dTelesco & Gezaril 119921: [Engelbracht et alj|2006l) . The surface 
brightness dynamic range between the starburst and the halo 
emission spans nearly three decades. Therefore, to subtract the 
central sources and obtain maps of the superwind and halo as 
minimally contaminated by the diffraction pattern of the bright 
center as possib le, we adapte d the CLEAN algorithm used in ra- 
dio astronomy (Hogbom ll974l) to substitute "clean" beams for 
beams with extended sidelobes. At each iteration, the peak of 
the map is found, and the scaled PRF at the same location is 
multiplied by a small gain and subtracted from the map; a nar- 
rower Gaussian with the same flux is added to a model map. The 
residual maps obtained after "cleaning" are shown in Fig.[TJ 

The starburst emission that we recover from the beam 
"cleaning" process, after adding back to the model maps the 
underlying residuals (that account for about 7% of the total), 
amounts to 337, 111, and 35.4 Jy at 250, 350, and 500 /mi, re- 
spectively, i.e. between 70% and 75% of the global emission. If 
we instead sum the emission from the four brightest sources se- 
lected from the 350 /mi SCUBA map, fitted with fixed positions 
and widths (~ 150pc) in all maps, then we already account for 
about 57% of the global emission in each SPIRE band. 

Since there is a smooth transition from the disk to the off- 
disk emission, the spatial separation between the two is some- 
what arbitrary. We note that the disk of M 82 is very thick and 
turbulent, and that the residual maps are still contaminated by 
disk emission. We thus masked pixels where the original flux 
density exceeds the cleaning threshold, as shown in Fig.fT] 

Dust opacities and emissivities: In all the following, we adopt 
the Galactic emissivity law derived by Reach et al. ( 1995) from 
the COBE survey. Dust grain opacities vary as a function of envi- 
ronment and are very uncertain. However, dust masses can eas- 
ily be rescaled if another value is assumed. We thus adopt an 
opacity valid for dust in the diffuse interstellar medium, k — 
1.62 m 2 kg" 1 in the DIRBE 140 /mi band (JLi & Draine 2001). 
Converting it to an opacity at 500 /mi, we obtain 0. 18 m 2 kg -1 , in 
good agreement wi th what would ha ve been obtained from the 
850 /mi opacity of lJames et all (120021) . i.e. (0.20 + 0.06) m 2 kg" 1 . 



3. The wind and tidal streams 

Cool dust is detected at 5cr significance in filaments out to 9 kpc 
from the center, at 250 /mi. The maximal extensions are to the 
north-east and to the south-east of the center, i.e. in filaments 
making a large angle with the polar axes of the galaxy and out- 
flow. As already remarked by lEngelbracht et aD (120061) for the 
PAH emission, dust occupies a very extended halo, and is not 
confined to the edges of what is tradit ionally cal l ed the super- 
wind. Comparison with the HI map of lYun et aD (1 19941) shows 
that the densest parts of the tidal streams are associated with 
some of the extended dust detected in SPIRE maps. Below, we 
provide quantitative arguments to support the impression that 
most of the dust mass outside the disk of M 82 has been removed 
by tidal interaction, rather than entrained by the superwind. 
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To approximately match the angular resolution of the maps 
and examine spatial variations in dust colors of the wind and 
halo, we convolved the "cleaned" maps with Gaussian kernels 
of appropriate widths, and regridded the resulting maps to the 
pixel size at the longest wavelength. Figure[2] shows the 250 ^/m 
to 350/im ratio map and the 250/im to 500 //m ratio map. Based 
on these ratios, the dust temperature in the wind and halo ranges 
between about 12 K and 50 K, and its average value decreases 
with projected distance from the starburst. Some filaments that 
are much warmer than their surroundings can be seen; they are 
located in-between bright dust spurs (Fig.|2]i, and thus must have 
relatively low densities. The present data do not allow us to test 
whether these color variations are caused purely by different ex- 
citation conditions, or whether they are affected by changes in 
the dust composition within the wind. Assuming that these fila- 
ments follow the Galactic emissivity law, we estimate that they 
have temperatures between 30 K and 50 K, i.e. at least as high 
as in the starburst. They seem to be associated with the starburst 
wind, rather than with the tidal features traced by the HI map. 

We estimate the total dust mass in the wind and halo from the 
500 fim and dust temperature maps, the latter derived from the 

250yum to 500yum ratio. We obtain M d halo = (l.l±0.5)xl0 6 M Q 
(mass error map summed linearly, and 2.5<x cut applied to the 
brightness maps). The mass derived from the global halo fluxes 
and a single-temperature fit (with Td ust - 22.3 to 25.6 K) is sim- 
ilar: M A h a i = (1.0 to 1 .2) x 10 6 M fT). This is much lower than 
the SCUBA estimate of I Alton et al] (U999I) (2 to 10 x 1O 6 M ), 
even using the same dust opacity. We cannot rule out a cold 
component undetected by SPIRE, but the discrepancy may pre- 
dominantly arise from several analysis issues: although the halo 
is barely detected in the SCUBA maps (because they are only 
2.6 kpc in diameter, and the chop throw placed the sky reference 
within the wind), the starburst and disk have not been subtracted 
with their sidelobes, and the halo emission may thus be severely 
contam inated by them; the procedure followed by lAlton et all 
( 1999) to estimate the dust temperature also entails large errors. 

The combination of the dust temperature and brightness 
maps shows that dust that is spatially coincident with the super- 
wind (likely to have been entrained by it) is warmer and more 
diffuse than dust that is associated with the tidal streams, the lat- 
ter consequently dominant in terms of mass, as illustrated by the 
dust mass map in Fig. [2] Defining large sectors to exclude the 
superwind regions, based on Ha and X-ray maps (with opening 
angles of 74° in the north and 64° in the south), we argue that at 
a minimum 65% of the halo dust mass is not associated with the 
superwind. Since tidal features overlap in projection with super- 
wind features, this fraction is a strict lower limit. 

The global dus t mass of M82, folding in the 1.25 mm mea- 
surement of lThuma et alJ (l2000h . is M d tot ~ (4.5 + 1 .4) x 10 6 M© 
(with Tdust = 24 K, which provides a fit too low by only 3% at 
350^m and by 22% at 250/rni). If our mass estimates are cor- 
rect, they imply that ~ 25% of all the dust has been expelled 
from the disk. Either the efficiency of gas and dust ejection by 
the superwind and the tidal interaction is very high, or the dust 
grain composition and size distribution in the wind and halo sig- 
nificantly differ from those in the starburst. Without data of com- 
parable quality at both shorter and longer wavelengths, we can- 
not test the second alternative, but we deem it less likely. 

Since the interstellar medium of M 82 must have become ex- 
tremely turbulent and porous as a result of the sustained injec- 
tion of mechanical energy by the successive starbursts, probably 
enough ultraviolet radiation is leaking out of the disk to provide 
most of the dust heating necessary to account for the observed 
emission. GALEX images show diffuse far-UV and near-UV 



light scattered by dust grains out to large distances dHoopes et alJ 
2005), at least up to 6 kpc on the southern side. However, this 
light is mainly confined to lobes around the minor axis, whereas 
dust emission is not. Optical photons must dominate the heating 
further away from the minor axis. 

Usin g the archival Spitzer 8 fim image to trace PAHs in the 
wind (Engelbracht et al. 2006), we find that the PAH emission 
does not extend as far out as the 250-500 fan emission, and that 
there is a steep radial gradient in the 8/im to 250 /im flux ratio 
(Fig. |2]l. There are also some regions of enhanced PAH bright- 
ness relatively close to the polar axis, which correspond well to 
the 250 //m spurs coincident with Ha emission. These regions 
contain relatively warm dust. We hypothesize that they are well 
outsid e the wind shock, wher e PAHs would have very short life- 
times (Micelot ta et al.l 12010). and that they simply correspond 
to more efficient illumination by UV photons leaking through 
discrete cavities in the disk. A large fraction of optical line emis- 
sion in t his region is also due to photoioniza tion, and not only to 
shocks dShopbell & Bland-Hawthornlll998"l) . The general gradi- 
ent in the 8 /mi to 250 //m flux ratio and in the far-infrared flux 
ratios would then be a natural consequence of the dilution of the 
radiation field with distance from the emitting stars. 



4. By way of conclusion 

The intergalactic medium can be enriched in metals by both tidal 
stripping and starburst winds. M 82 offers an exquisite example 
of both processes at work in the same galaxy. Considering the 
spatial distributions of the dust brightness and temperature in 
the halo, the tidal interaction may be far more effective than the 
wind in removing the cool interstellar medium from the disk of 
M 82. We estimate that more than two thirds of the extraplanar 
dust mass is associated with the tidal streams traced by atomic 
hydrogen. If dust is able to survive for a long time in the ex- 
tragalactic medium around such interacting systems, which may 
have been common in the past, then it will have significant col- 
umn densities along lines of sight to distant galaxies, and cause 
extra attenuation outside the Milky Way. 
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Fig. 1. SPIRE maps resulting from the "clean ing" process. Top row: residual maps showing the wind an d halo emission. Contours 
drawn from the HI map of I Yun et al.l (1994) and from the Ha- [Nil] map of Boselli & Gavazzi (2002) are superimposed on the 
500 /im and 350jum maps, respectively. Bottom row: subtracted central emission. The stellar disk is delineated by the IRAC 3.6//m 
contours on the 500 /im map. The initial maps have peak brightnesses of 105, 48.2, and 21.6 Jy/beam, at 250, 350, and 500 jum 
respectively. The beam FWHM is sketched to the top right. 
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Fig. 2. Top row: Flux ratio maps. Left: 250//m to 350/im ratio, at 28" angular resolution (FWHM). Right: 250/mi to 500 fim ratio, 
at 40" angular resolution (FWHM). Contours of the 250/mi map are overlaid in black. Bottom row: Left: dust mass map in Mq per 
pixel (FWHM of 40", corresponding to 760 pc; pixels are 170pc on a side). HI contours are superimposed in white and Ha contours 
in black. Special note: the dust blob located around 09:55:23 +69:46:10 is likely associated with M82, rather than a background 
source, since its morphology in all SPIRE bands is diffuse. Right: 8//m to 250/mi flux ratio map of the wind and halo, at the 250//m 



